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Over  the  past  two  decades  passive  microwave 
imaging  systems  have  proven  to  be  effective  recon¬ 
naissance  tools  in  polar  environments.  However, 
the  mechanical  scan  mechanism  and  high  gain 
electronics  characteristic  of  this  class  of  sensors 
commonly  impart  noise  and  unwanted  artifacts 
to  image  data  they  produce,  complicating  visual 
analysis  and  automated  classification  procedures. 
The  fact  that  data  in  individual  scan  lines  are 
characterized  by  statistical  stationarity  and  that 
information  in  adjacent  pixels  is  highly  correlated 
due  to  oversampling  of  these  image  data  suggests 
that  Wiener  multichannel  filtering  techniques  may 
prove  effective  in  this  application.  Wiener  filters 
applied  to  pa.ssicc  microwave  images  of  first-year 
.sea  ice  were  constructed.  Four  mopr  parameters 
that  define  the  filter  (lag  or  pixel  offset  between  the 
original  and  desired  scenes,  filter  length,  number 
of  l.  xes  in  the  filter,  and  weight  applied  to  the 
empirical  correlation  functions)  were  varied.  Re- 
sults  were  compared  visually  to  asse.ss  the  effect  of 
each  variable  on  image  quality.  Effective  filters 
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that  limit  high  frecfuency  noise  and  enhance  ice 
characteristics  use  a  lag  of  one  pixel,  consist  of  two 
or  three  channels,  are  five  pixels  in  length,  and 
weigfi  the  auto-  and  cross-correlation  functions 
equally. 


INTRODUCTION 

Passive  microwave  imaging  systems  represent  one 
component  of  a  balanced,  multisensor  instrument 
suite  that  is  required  to  characterize  the  sea  ice 
cover  in  polar  regions.  A  series  of  field  e.^peri- 
ments  extending  over  almost  2  decades  involving 
surface  investigations  and  aircraft  overflights  es¬ 
tablished  the  utility  of  passive  microwave  data  in 
assessing  the  relative  age  of  sea  ice,  its  geographic 
extent  and  concentration,  the  location  of  water 
openings,  and  the  dynamic  structure  of  the  polar 
ice  pack  (Campbell  et  al.,  1975;  1976;  1978;  Hol- 
linger  et  al.,  1984;  Eppler  et  al.,  1986;  Cavalieri 
et  al.,  1986).  The  overall  results  of  these  missions 
demonstrate  that  aircraft  passive  microwave  im¬ 
ages  of  high  cjuaiity  are  of  significant  value  to 
programs  that  involve  reconnaissance,  operations, 
and  research  in  polar  environments. 
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Aircraft  sensors  constructed  to  provide  fine 
resolution  passive  microwave  inuiserv-  typically 
employ  a  scanning  antenna  assembly  to  measure 
the  extremely  weak  (10  VV)  microwave  signal 
emanating  from  the  earth’s  surface  (Hollinger  et 
al.,  1976;  Eppler  and  Heydlauff,  1990).  High  gain 


electronics  amplify  the  signal  and  pass  it  to  re¬ 
cording  and  display  systems  in  the  aircraft  eahin. 
Both  the  mechanical  scan  mechanism  and  high 
gain  electronics  commonly  impart  imdesirahle  ar¬ 
tifacts  to  the  data  that  complicate  Ijoth  automated 
classification  procedures  applied  to  derive  ice 


Figure  1.  KRM.S  image  of  first-year  sea  ice.  The  scene  is  .512  x. 512  pixels  in  size  anti  shows  an  area  approxiniatcK  1-1.5 
km  across.  The  image  was  aciiuirecl  from  an  altitmle  ol  6100  in  (20.000  ft)  on  1 1  March  1988  in  the  B<-aufort  Sea  200  kin 
from  the  Alaskan  coast  near  71°45'N  144°30'\\5  Surfaces  shown  range  in  Isrightness  temperature  from  approximate))  215 
K  to  242  K.  Dark  areas  are  radiomefricall)'  warm  and  indicate  ice  that  i.i  highly  emissive  at  the  .sensor  Iretiuency  (3.3  6 
GHz):  light  areas  are  radiometrically  cool.  The  mottled  pattern  is  Ixpical  of  hrst-year  ice  in  this  region  and  rellects  varia¬ 
tion  in  surface  properties  of  the  ice  surface  and  overlving  snow. 
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type  distributions  (Holyer,  1990;  Eppler  et  al., 
1986;  Eppler  and  Farmer,  1991b)  and,  to  a  lesser 
degree,  visual  analysis  of  the  images.  Electronic 
artifacts,  which  arise  from  noise  introduced  when 
the  weak  radiometer  signal  is  amplified,  impart  a 
snowy  appearance  to  the  image  (Fig.  1).  This  adds 
to  variance  associated  with  intensities  typical  of 
different  surface  types  and  increases  the  probabil¬ 
ity  that  pi.xels  will  be  misclassified  when  segmen¬ 
tation  techniques  are  applied  (Eppler  et  al.,  1986), 
especially  those  associated  with  automated  classi¬ 
fication  systems.  Scan  artifacts,  which  consist  of 
magnitude  shifts  in  signal  intensity  and  in  lateral 
shifts  in  feature  location  from  scan  to  scan,  arise 
from  imbalance  betw’een  antennas  and  from  tim¬ 
ing  aberrations  that  occur  during  data  transfer 
and  storage.  This  results  in  each  scan  being  classi¬ 
fied  as  an  edge  when  texture  classifiers  based  on 
gradient  and  edge  detection  logic  are  applied  to 
the  image  (Eppler  and  Farmer,  1991b). 

Work  described  here  reports  results  of  apply¬ 
ing  multichannel  Wiener  filtering  techniques  to 
minimize  electronic  noise  imparted  to  fine  resolu¬ 
tion  passive  microwave  sea  ice  images  by  high- 
gain  sensor  electronics.  In  particular,  we  apply 
the  method  to  images  of  first-year  sea  ice,  which 
show  subtle  patterns  of  variation  in  surface  emis- 
sivity  (Figs.  1  and  2).  These  patterns,  which  are 
particularly  susceptible  to  obscuration  by  instru¬ 


ment  noise,  are  important  from  a  scientific  per¬ 
spective  because  they  reflect  areal  variation  in 
physical  characteristics  of  snow  and  ice. 

Image  data  used  to  evaluate  the  filter  were 
obtained  with  the  U.S.  Navy  K^-band  Radiometric 
Mapping  System  (KRMS),  an  aircraft  sensor  that 
acquires  high  resolution  passive  microwave  imag¬ 
ery  at  a  frequency  of  33.6  GHz  (Eppler  et  al., 
1984;  1986;  Eppler  and  Heydlauff,  1990).  Our 
objectives  in  applying  the  filter  are  to  improve 
visual  clarity  of  the  images  so  that  ice  variability 
may  be  mapped,  and  to  enhance  the  effecfiveness 
of  automated  image  classification  techni({ues  by 
noise  reduction.  As  will  be  shown  belovs’,  use  of 
Wiener  signal  enhancement  filters  is  appropriate 
in  this  application  as  a  natural  consequence  of  the 
structure  of  data  acquired  with  a  scanning  sensor. 

BACKGROUND 

Traditional  methods  of  dealing  with  noise  reduc¬ 
tion  use  various  smoothing  filters  oi  applications 
of  the  Fourier  transform  functions.  Several  of  the 
traditional  techniques  were  applied  to  passive 
microwave  imagery,  largely  with  unsatisfactory 
results.  For  example,  consider  the  results  of 
applying  a  moving  average  function  (convolution) 
to  the  three  images  of  first-year  sea  ice  in  Figure 


Figure  2.  KRMS  test  images  used  in  this  study.  Two  of  the  images  (B  and  C)  show  Krst- 
year  sea  ice  exclusively.  The  third  image  ^A)  consists  predominantly  of  first-year  ice.  hut 
also  shows  a  break  in  the  ice  cover  (bright  streak  at  center  of  image  in  first  column)  that 
exposes  sea  water,  which  is  radiomcti  ically  cold.  This  image  (A)  was  included  to  test  the 
degree  to  which  the  Wiener  filter  technique  is  robust  with  respect  to  high-contrast  fea¬ 
tures  that  occur  in  an  otherwise  low-contrast  .scene.  Each  image  is  128x  128  pixels  in 
size.  The  test  images  are  composed  of  scans  acquired  w’ith  only  one  of  KRMS's  three  an¬ 
tennas  and  so  are  foreshortened  in  the  downtrack  direction.  The  area  shown  b>-  each  test 
image  is  approximately  3.6  km  across  and  10.9  km  from  top  to  Inittom.  Histogram  eipial- 
ization  and  range  transformation  were  also  performed  over  all  three  images. 
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2.  Three  moving  average  functions,  2x2,  3  x  ,3, 
and  5x5,  were  used  to  produce  the  filtered  im¬ 
ages  in  Figure  3.  Convolution  reduces  noise  and 
eliminates  the  snowy  appearance  of  the  original 
image,  but  at  the  expense  of  high  frequency  infor¬ 
mation.  Information  removed  by  the  filter  can  be 
displayed  in  image  torm  by  subtracting  the  filtered 
image  from  the  original  image  (Fig.  4).  These 
image  difference  maps  are  convenient  for  visually 
presenting  the  net  change  in  an  image  due  to  a 
filter  operator  (original  image  =  coinolved  image 
+  difference  image).  Note  in  Figure  3  that  fea¬ 
tures  become  fuzzy  as  a  result  of  applying  these 


filters.  This  blurring  is  unacceptable  la  t  ause  tex¬ 
tural  and  edge  information  is  being  reinox  ed  (Fig. 
4)  that  might  be  useful  in  subse<juent  analytical 
processing,  f'igure  4  captures  these  regions  of 
greatest  smoothing  as  ghostlike  features  on  an 
otherwise  featureless  background.  W'v  seek  to 
preseiwe  this  high-frequency  information  while 
removing  the  noise,  and  standard  convolution 
methods  (Rush,  1990)  have  failed  with  this  data, 
.Application  of  other  filters  such  as  the  double 
Fourier  transfo*-m  (Gonzales  aiul  \\  intz.  1977; 
Rush,  1990)  failed  also  to  proxide  the  desired 
enhancc'inent,  probably  due  to  the  fact  that  noise 


Figure  3.  Images  in  Figure  2  (labele<l  original  in  tl.is  aiul  siibseqiieiil  hgtires)  eoiivoKed 
with  2x2,  ,3  X  .3,  and  .5  x  .5  smoothing  kernels.  Note  the  increased  hhirring  ot  subtle  lev- 
tural  features  and  the  obscuration  of  feature  edges  by  tiperators  of  increasing  si/e.  .Ml  tbi 
images  in  Figures  3  and  4,  along  with  all  subsequent  images,  are  presented  bistograin 
equalized  and  range  e.xpanded  to  maximize  the  visual  impact  of  subtle  features. 
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present  in  the  images  is  not  arrayed  within  dis¬ 
tinct  spatial  frequencies  which  are  necessary  to 
define  elfective  Fourier  filters.  Other  less  popular 
approaches  were  tried  with  similarly  disappoint¬ 
ing  lack  of  success. 

\\1ENER  ENHANCEMENT  FILTERING 

iMost  current  theories  for  prediction  and  enhance¬ 
ment  of  time  series  signals  are  based  on  the  origi¬ 
nal  work  of  Norbert  Wiener  (Wiener,  1942).  Wie¬ 
ner  was  interested  in  the  systematic  design  of 
filters  to  enhance  signals  and  minimize  unwanted 
noise  in  time  series  analysis.  In  particular,  Wiener 
desired  to  design  filters  that  could  be  used  for 
signal  prediction  and  forecasting.  Some  of  the 
earliest  applications  of  Wiener’s  work  involved 
long-range  weather  forecasting  during  the  Second 
World  War.  Most  current  applications  of  this  form 
of  analysis  concern  separation  of  desired  signals 
from  unwanted  noise  in  seismic  signal  analysis 
(Robinson,  1983). 

Examples  of  the  extension  of  Wiener  theory  to 
multiple  channels  are  numerous  and  span  diverse 
research  fields.  One  of  the  more  recent  nonseis- 
mic  applications  is  reported  in  Williams  et  al. 
(1988),  who  discuss  applications  for  multivariate 
data  including  atmospheric  isotopic  variations 
through  earth  history.  In  multivariate  applications 
such  as  these,  each  of  the  channels  is  assumed  to 
record  similar  events  in  time  which  have  been 
corrupted  by  noise  that  is  uncorrelatable  across 
channels.  Wiener  filters  are  applied  to  enhance 
correlative  events  and  subdue  noise  manifest  by 
uncorrelatable  events  observed  across  channels. 
Wiener-type  theory  also  has  been  widely  used  for 
image  restoration  (Gonzalez  and  Wintz,  1977; 
Rosenfeld  and  Kak,  1982;  Rush,  1990),  although 
the  manner  in  which  we  apply  these  filters  here 
is  derived  more  from  seismic  literature  (Robinson 
and  Treitel,  1980;  Robinson,  1983;  Robinson  and 
Durrani,  1986)  and  its  extension  into  the  multi¬ 
channel  domain  (Williams  et  al.,  1988)  than  from 
conventional  methods  of  image  restoration. 

In  the  current  application  we  consider  imag¬ 
ery  acquired  by  scanning  sensors  to  be  analogous 
to  time  series  data.  Rows  or  columns  of  pixel 
intensities,  which  represent  time-dependent  mea¬ 
surements  acquired  by  the  sensor  antenna  as  it 
sweeps  across  the  imaged  field  of  view,  form  the 


time-varying  signal  to  be  filtered.  .Adjacent  scans 
(rows  or  columns)  contain  similar  information  that 
is  highly  correlated  for  two  reasons.  First,  adja¬ 
cent  scans  carry  redundant  information  because 
sciumer  data,  if  acquired  properly,  are  o\  en>ainpled 
to  satisfy  Nyquist  criteria.  Second,  different  parts 
of  features  that  span  more  than  one  scan  appear 
at  similar  points  in  adjacent  scans  (Fig.  5).  Appli¬ 
cation  of  Wiener  multichannel  techni(iues  pro¬ 
vides  a  means  of  a)  enhancing  edges,  features, 
and  patterns  that  appear  in  adjacent  scans  and 
that  are  highly  correlated  with  respect  to  position 
in  the  scan  and  b)  minimizing  the  effect  of  noise 
that  is  not  correlated  between  scans. 

Multichannel  Wiener  filtering  is  described 
elsewhere  in  detail  (Rosenfeld  and  Kak,  1982; 
Robinson,  1983).  Principles  that  underlie  this  type 
of  filtering  are  summarized  below'  so  that  assump¬ 
tions  inherent  in  correct  application  of  the  method 
are  clear.  The  notation  of  Robinson  (1983)  will 
be  used  in  this  discussion,  and  it  will  be  cast  in 
terms  of  image  data  addressed  by  this  study. 

Wiener  filtering  procedures  are  based  on 
least-squares  techniques.  The  least-squares  class 
of  filters  depends  on  the  existence  of  correlations 
between  the  actual  and  the  desired  scene  to  de¬ 
fine  an  optimal  filter,  which  is  used  to  convoKe 
the  actual  scene  into  the  desired  scene.  Let  us 
define  the  input  scene  as  X^n  and  the  desired 
output  as  Z„i,„  where  m  represents  the  time- 
sampled  raster  scan  (row'  or  column  of  pixels)  and 
n  represents  the  number  of  scan  lines  that  w'e  wish 
to  be  included  in  the  Wiener  filter.  Correlations 
should  exist  among  the  information  expressed  by 
the  family  of  n  scan-lines  chosen.  Note  that  if  the 
noise  spectra  w'ere  correlated  with  the  signal, 
then  that  portion  of  the  noise  would  be  considered 
as  signal  in  this  analysis,  potentially  amplifying 
undesired  effects.  However,  as  noted  in  the  next 
section,  the  noise  we  are  trying  to  remove  is  not 
correlatable  to  the  signal  in  the  images  used  in 
this  study.  We  define  the  standard  autocorrelation 
function  of  the  input  for  each  scan  line  n  as 

4>„/T)  =  £{x„,jX’;„-,)J,  r  =  0,l,2,...,  (1) 

where  (r)  represents  the  autocorrelation  func¬ 
tion  of  raster  scan  J  (I  n)  over  m  pixels  and 
for  pixel  lags  (t)  0,1,2,...  and  so  forth,  and  E 
represents  the  usual  recipe  for  calculating  the 
autocorrelation  of  raster  scan  x„y  in  terms  of  the 
line  itself  and  its  transpose  defined  at  various  lags. 
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The  T  represents  the  standard  transptvse  function. 
In  our  case,  m  4-  t  =  total  number  of  pixels  in  the 
original  raster  scan.  That  is,  X„.„  represents  the 
image  defined  by  pixels  1,  2,  ...,  m  for  the  first  n 
scan  lines  while  Z„,„  is  the  image  defined  by  pixels, 
T,  r+1,  ...,  m  +  T,  for  the  .same  n  .scan  lines  of 
the  original  image.  Similarly,  the  cross-correlation 
function  can  be  defined  as 


=  r  =  0, 1,2,3,...,  (2) 


where  ~  refers  to  the  corresponding  raster  scan 
in  the  desired  image. 

These  autocorrelation  and  cross-correlation 
functions  form  the  basis  from  w'hich  a  filter  is 
derived  that  best,  in  a  least-squares  sense,  con¬ 
volves  the  input  scene  into  the  output  scene.  That 
is,  if  the  best-fit  filter  for  raster  scan  k  (where 
1  <  <  n)  is  defined  as 

fkn{s)fku{s)  ■  ■  ■ 

•  •  • 

/kn\{s)fk2lis)  ■  ■  ‘/w(i')J 

where  m  and  n  are  defined  as  before,  p  is  the 
length  of  the  filter,  and  s  =  0,  1,  2,  ....  (p  -  1),  then 
each  element  tjij  of  the  actual  output  scene  tj  is 
given  by 


/iv=/i(s)  = 


(4) 

where  i  =  1,  2,  ...,  m,  j=  1,  2,  ...,  n,  k  =  I,  2,  ...,  n. 
The  error  between  the  desired  output  and  the 
actual  output  is  simply  the  difference  between 
the  desired  output  and  the  actual  output.  It  is  this 
error,  normalized  in  the  usual  sense,  that  we  wish 
to  minimize  in  a  least-squares  algorithm.  Each  set 
of  filter  coefficients  (ft,)  can  be  obtained  by  solving 
p  sets  of  simultaneous  equations  given  by 

/*«<I)„*(0)  +/*,<£>„*(  -  1)  -!  •  •  •  +/*,<*>„<(  -  s)  =  <D,40). 

+  fki^uki^)  +  •  •  •  +  fh<Py,k{  1  -  .5)  = 

•  •  • 

fw^xik{s)  +  -  1 )  +  ■  •  •  +fk,^xyk{0)  =  ‘I>i,A(.s),  (.5) 

where  1  <  A:  <  n.  These  sets  of  normal  equations 
form  a  block  Toplitz  matrix  and  can  be  solved  with 
relative  ease  using  one  of  a  number  of  efficient 
algorithms.  Robinson  (1983)  provides  an  excellent 
discussion  of  one  such  method. 


APPUCABILITY  OF  WIENER  THEORY  TO 
IMAGE  DATA 

Applicability  of  multichannel  Wiener  filtering  to 
the  current  data  set  is  dependent  upon  a  series  of 


assumptions  that  must  l)e  shown  to  be  reasonable 
with  respect  to  image  data.  .Xs  Rof)inson  (19831 
notes,  the  extent  to  which  results  of  the  filter  can 
be  applied  depends  upon  the  <  xtcnt  to  which 
these  assumptions  are  satisfied  by  the  input  data. 
Four  assumptions  are  considered  here:  1)  that 
the  image  data  analyzed  cotdorjn  to  statistical 
stationarity  (that  is,  statistical  properties  of  the 
data  do  not  vary  greatly  over  the  scan),  2)  that 
mean  square  error  between  the  input  and  output 
sequences  represents  an  optimal  criterion  to  eval¬ 
uate  the  filter,  3)  that  the  derived  filter  is  a  linear 
operation  and  is  of  finite  length,  and  4)  that  adja¬ 
cent  pi.xels  in  the  image  in  fact  displav  some 
degree  of  statistical  dependence.  Before  these 
assumptions  are  discussed,  however,  the  nature 
of  the  noise  associated  with  these  images  needs 
to  be  addressed.  The  two  major  .sources  of  noise 
are  the  noise  associated  with  the  amplifier  and 
the  thermal  noise  associated  with  the  surface  sig¬ 
nal.  Both  of  these  sources  of  noise  are  uncorre¬ 
lated  with  signal  and  displays  purely  additive  char¬ 
acteristics  (as  opposed  to  multiplicative).  For 
example,  examination  of  scans  that  cross  low  in¬ 
tensity  and  high  intensity  show  similar  variance 
structure  when  e.xamined  visually.  If  noi.se  were 
multiplicative,  the  variance  structure  would  be 
expected  to  be  dissimilar  over  these  features. 

Assumption  1 

The  first  assumption  is  that  the  statistical  proper¬ 
ties  of  the  input  sequence  X„,„  and  the  desired 
output  sequence  Z,„„  (which  in  this  case  consist 
of  rows  and  columns  of  pi.xels  in  the  unfiltered 
and  the  filtered  passive  microwave  scene)  do  not 
change  with  position  within  the  portion  of  the 
scene  being  analyzed.  This  assumption,  which  is 
an  adaptation  of  the  theory  of  statistical  stationar¬ 
ity  applied  commonly  to  time  series,  implies  that 
the  probability  of  encountering  any  specific  inten¬ 
sity  value  at  a  pixel  location  is  not  tied  to  a 
specific  origin  and  that  the  scene  being  analyzed 
is  conceived  to  be  effectively  infinite  in  extent 
(adapted  from  Wadsworth  et  al.,  1953).  Visual 
examination  of  individual  raster  scans  of  the  input 
test  images  (Figs.  2  and  5)  presented  in  this  study 
suggests  that  the  range  of  variability  is  reasonably 
uniform  both  along  individual  rows  (along  image 
scans)  and  dow'ii  successive  image  columns.  Statis¬ 
tics  that  describe  variation  in  pixel  intensity 
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Figure  5.  Plots  of  scans  300-309.  each  512  pixels  wide,  troin  the  KKMS  iinaKe  shown  in 
Figure  1,  Peaks  and  valleys  that  represent  major  features  can  he  followed  readily  iroin 
scan  to  scan. 


within  and  between  scans  support  this  conclusion 
(Table.s  1  and  2). 

Assumption  2 

The  second  assumption  is  that  the  mean  .stpiare 
error  between  A,„„  and  Z,„„  is  the  optimal  criterion 
to  evaluate  performance  of  the  filter  operators 
determined  in  the  analysis.  This  is  similar 
to  assumptions  associated  with  regre.ssion  analy¬ 
sis,  principal  component  analysis,  and  lea.st- 


squares  procedures  in  general,  (’onditions  under 
which  tliis  assumption  is  violated  occur  either 
when  the  statistical  variability  is  correlated  spa¬ 
tially  within  an  image,  or  when  the  autocorrela¬ 
tion  function  e.xhibits  a  paucity  of  correlatable 
events.  Given  the  charactt'r  of  image  rlata  eonsirl- 
ered  here,  we  find  no  ei)inpelling  reason  why 
this  assumption  .should  not  hold.  Tables  1  and  2 
inrlicate  that  the  data  are  not  spatially  correlated. 
Additionally,  Nyquist  sampling  used  to  acajuire 
the  data  produces  overlap  in  information  con- 
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Table  1.  Statistics  lor  Radiances  in  Raiidointy  Selected 
Image  ('oliimns 


Columti 

Miniittuin 

Radiance 

fKj 

Maximum 

Radiance 

(E) 

Mean 

Radiance 

(Ej 

Stamlard 

Deviation 

9 

220.0 

2;55.0 

22.S.  1 

2.0 

13 

221,1 

233..S 

22.S.4 

2.1 

39 

219.3 

233.5 

220.5 

2  S 

43 

221.7 

233.0 

220.7 

2,3 

HO 

220. 1 

2.34.1 

220.0 

2-4 

f>2 

221.2 

233.2 

2202 

23 

(iS 

219.5 

231.1 

220.7 

2.5 

74 

210.4 

232.0 

220.2 

2  1 

SI 

219.2 

2.30.0 

225.9 

2  1 

S2 

210.9 

230.0 

225  7 

2  3 

H9 

214, .5 

232.3 

225.5 

2.5 

99 

21.5.0 

230.3 

225,0 

2  7 

101 

217  4 

230.9 

224..S 

2.0 

109 

210.9 

231.7 

224.3 

2  S 

120 

219.1 

231  1 

224  5 

2  0 

tained  in  adjact?nt  pixels  (Fig.  6lj)  (Eppler  and 
Heydlauff,  1990)  and  ensures  that  events  re¬ 
corded  in  adjacent  scans  will  be  highly  correlated. 

Assumption  3 

The  third  assumption  is  that  the  optimal  filter 
derived  by  the  least-scjuare  analysis  is  a  linear 
operation  and  is  of  finite  length.  The  validity 
of  this  assumption  relative  to  passive  microwave 
measurements  cannot  be  completely  demon¬ 
strated.  This  assumption  probably  is  valid  for  the 
relatively  simple  images  of  first-year  sea  ice  to 


Table  2.  .Statistics  tor  Railianccs 
Image  Rows 

ill  Raiidonily 

Selecteil 

Row 

Minimum 

Radiance 

(E) 

Maximum 

Radiance 

(E) 

Mean 

Radiance 

(E) 

Standard 

Deviation 

3 

219.5 

233.2 

220.4 

2.9 

13 

219.3 

231.4 

227.5 

2.5 

40 

222.7 

235.0 

229.0 

2  4 

4<) 

219.0 

233.7 

227  1 

3.3 

50 

220.3 

231.2 

220.9 

2.5 

70 

220.7 

232.2 

220.1 

2.5 

.S.3 

221.5 

2.30..S 

220.1 

2.2 

91 

220.0 

231.7 

225.0 

2  1 

92 

221  5 

2.30.0 

■220.2 

2.0 

100 

222.5 

232.3 

220.5 

2.4 

103 

220.7 

229.2 

225.3 

2.1 

107 

220.2 

229.7 

225.1 

1.7 

lO.S 

220.5 

2.30.H 

22.5.4 

2.1 

1  15 

220.5 

230.5 

225.5 

2.2 

124 

215.S 

2.32.9 

225.4 

.3.4 

which  we  apply  the  method  here,  (.'oncr'i  abic, 
nonlint'ar  operators  may  prose  to  be  of  greater 
usefulness  wlum  complex  scenes  that  iticlude 
multiple  surface  types  are  anal\/,ed,  esptcialK 
those  scenes  in  which  larger-scale  \ariabilit\  oxer 
short  distances  is  juxtaposed  with  (and  indepen¬ 
dent  of)  small-scale  textural  patterns. 

Assumption  4 

The  final  assumption,  and  probably  tin*  most  im¬ 
portant  one,  is  that  adjacent  pixels  display  some 
degree  of  statistical  {lependenc<‘.  l  )ei)endence  })e- 
tween  successive*  obserx  ations  is  one  of  the  funda¬ 
mental  characteristic  of  time  si'rics  data  (Wads¬ 
worth  et  ah,  1953).  .Application  of  Wiener 
techniques  to  data  that  lack  such  correlation  pro¬ 
duces  unpredictable  results.  The  inherent  struc¬ 
ture  and  scale  of  image  data  anaK  zed  here  guaran¬ 
tees  that  a  certain  degree  of  dependence  exists 
between  adjacent  pixels,  both  b(*cau.sc  the  data 
are  sampled  robustly  with  respect  to  overlap  (Fp- 
pler  and  Heydlauff.  1990)  and  because  features 
that  define  the  signal  extend  oxer  areas  of  the 
image  large  enough  to  encompass  multiple  pixels. 

Each  observation  (pixel)  rc'presents  the  mea¬ 
surement  of  microxvaxe  emission  integrated  oxer 
•in  elliptical  area,  the  size  and  shape  of  which 
depends  on  viewing  angle  of  the  si*nsor  (Fig.  6a). 
\  significant  fraction  of  the  information  carried 
in  adjacent  pixels  is  redundant  because*  raxx  data 
are  collected  and  sampled  such  that  footprints  of 
adjacent  pixels  overlap  (Eppler  and  Heydlauff, 
1990).  The  radiometer  signal  is  oxersamplc  d  in 
the  alotig-scan  direction  (across  image  roxvs)  such 
that  the  footprint  of  each  pixel  overlaps  that  of  its 
neighbor  by  approximately  H()9(-  (Fig.  6b).  Ciround 
.swaths  of  adjacent  scans  in  test  images  used  in 
this  study  overlap  by  approAimately  40%,  mean¬ 
ing  that  the  footprint  of  adjacent  pixels  overlap 
hy  a  similar  amount  in  the  across-scan  direction 
(down  image  columns). 


METHODOLOGY 
KRMS  Characteristics 

Imagery  used  here  xxas  aciiuired  in  March  19(SS 
in  the  Beaufort  St'a  n'gion  of  coastal  .Alaska  using 
the  K,,-band  Hadionietric  .Mapping  .System  (KRMS) 
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Flight  Direction  — ^ 

Fistin’  li  A'j  Variation  in  the  size  and  shape  ol  tlie  KHMS  3  rlB  hearnspol  ,is  a  fnnc  tion  ol 
look  auj;le.  The  beam,  which  subtends  a  eirenlar  spot  at  nadir,  subtends  an  elliptie.d  spot 
at  off  nadir  angles.  Beamsiiots  are  plotted  to  scale.  Two  conversion  lactors  .ire  giscn  lor 
each  angle.  The  left  nninher.  when  nuiUiplied  by  aircra.l  altitude,  gives  the  area  of  the 
beamspot  in  sijnare  units  of  altitude.  The  area  subtended  by  the  nailir  bcainspot  from 
10, OOP  ft  altitude  thus  is  2302  sq  ft.  Tin-  right  inniihcr  r<-presents  the  relatirc  ,uca  ol  the 
beamspot  with  respect  to  the  area  of  the  nadir  In-aiiispot  The  .irea  imaged  by  c.ich  beam 
at  hO”  thus  is  3,73.39  tinn-s  larger  than  the  area  imaged  at  nadir  B)  Overlap  between 
beamspots  in  adjacent  KHM.S  scans  plotted  for  diHerr-nt  hxik  angles,  KKMS  antenna  scan 
rate  is  adjusted  with  respect  to  downtrack  aircraft  motion  as  rlata  are  l  ollected  such  that 
the  beam  in  each  successive  scan  'veeps  a  ground  swath  that,  at  nadir,  iiu  ludes  SD'^r  ol 
the  area  imag»-d  in  the  previous  scan.  Overlap  at  off  nadir  angles  is  grc.iler  than  SO'f  be¬ 
cause  the  size  of  the  beamspot  increases  with  look  angle. 
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in  conjunction  with  the  Na\  y-NASA  SSM  / 1  Vali¬ 
dation  Experiment  (Fanner  et  ah,  1989;  Eppler 
and  Farmer,  1991a;  Cavaheri  et  ah,  1991).  KlhMS 
is  an  airborne  passive  microwave  imaging  system 
that  senses  \ertica!ly  polarized  radiation  emanat¬ 
ing  from  the  Earth’s  surtace  ai  a  center  Irecjneney 
of  .33  6  GHz.  In  its  present  confignratiom  the 
KRMS  antenna  assembly  is  enclosed  in  a  pod  that 
hangs  either  from  the  bomb-bay  or  wing  of  a  F-3 
class  aircraft.  Three  parabolic  antennas,  mounted 
120°  apart  on  a  single  shaft  within  the  pod,  rotate 
about  a  horizontal  axis  that  is  parallel  to  the 
direction  of  flight.  Only  the  signal  from  the  down- 
ward-facing  antenna  is  sampled  and  recor<led  at 
any  given  time.  Scan  is  across  track  in  a  vi'rtic'a) 
plane  and  provides  coverage  of  .30°.  both  sides  of 
nadir  (Fig.  7).  The  rate  at  which  the  antenna 
assembly  rotates  is  controlled  by  the  operator  in 
the  aircraft  cabin  and  is  set  based  on  altitude  and 
ground  speed  to  provide  Hi)%  overlap  between 


successive  scans.  Antenna  i)eam  vvidtis  is  1°,  so 
ttiat  spatuil  resolution  of  the  unprocessed  sigmd 
is  approximately  16  ft  per  1000  ft  of  flight  altitude 
at  nadir.  Kadiometric  sensitivity  measured  in  the 
laboratory  is  (1.05  K/s;  operational  sensitivity  is 
estimatetl  to  Ire  0.5  K  or  better.  Additional  infor¬ 
mation  concerning  technical  characteristics  of  tlie 
KRM.S  sensor  is  provided  by  Eppler  et  al.  (1984; 
1986). 

Raw  data  are  ix'cordt'd  on  anal<«;  tape  and 
convt'rted  to  digital  form  aft<*r  the  aircraft  mission 
is  complete  (Eppler  and  Heydiauff,  1990).  Digi¬ 
tized  images  are  .512  pixels  wide.  Analog  data 
are  sampled  to  provide  80%  overlap  between 
adjacent  [rixels.  Each  row  of  512  pixels  corre¬ 
sponds  to  oiu'  scan  by  a  single  antenna  across  the 
full  100°  field  of  view.  Suceessive  lines  contain 
radiance's  derived  from  a  different  one  of  the  three 
antennas  which  repeat  in  setjuence  (i.e.,  0,  1,  2. 
0.  1 ,  2,  0,  1 ,  and  so  forth). 


Fitiuri'  7.  can  geometry  of  the  K.RM.S  sensor.  The  dovvnwarcl-looking  antenna  asseiniily 
rotates  and  scans  a  l()0°-\vide  swath  beneath  the  aircraft  centered  on  nadir.  Forward  mo¬ 
tion  of  the  aircraft  moves  the  footprint  of  successive  swaths  farther  along  track.  Rotation 
rate  of  the  antenna  is  set  such  that  one  swath  overlaps  the  next  by  SOr?-  at  nadir.  The  inci 
dcnce  angle  from  which  the  antenna  views  the  surface  varies  witli  position  along  a  scan. 
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Test  Data 

'rhrt-f  (fst  ittia.iifs  <»l  Hrst->far  sea  iee.  each 
128  X  128  pixels  in  si/.e  (Fi^.  2),  xsere  extraefeci 
fiorn  the  eential  portion  ol  a  full  KHMS  scene 
{512x  512  pixels).  These  smaller  128  x  128-pixel 
images  \\«‘re  usetl  instea<l  ol  the  hill  512-pixel- 
wide  iinasies,  both  to  rednee  computation  time 
and  to  minimize'  eross-traek  xarianee  in  miero- 
xxaxt'  radiance.  Mierowas f  emissixitx  \aries  as  a 
hinetion  ol  fhc'  incidc'nee  anulc'  irom  which  a 
surface'  is  xiexxc'd.  Fart  ol  the  xariancc'  ohserxt'd 
in  KRMS  data  results  Irom  cliani'c's  in  look  an^lc' 
from  pixc'l  to  pixc'l  across  the'  100“  held  of  xie'w. 
(aiolc'r  hrightiu'ss  te-mperaturc'S  (lighter  tones)  oh- 
sc'rxc'd  down  the  ee'iitral  section  of  Figure'  1  as 
eomparc'cl  with  warmer  temperatures  (darkc'r 
tone's)  at  image  limbs  re'Snlt  Irom  cliHerc'iiec's  in 
the  angle  at  whic  h  KRMS  antennas  view  the  ice 
surface'  (Fig,  7).  Theorx  predicts  and  hc'lcl  expc'ri- 
m  *nts  shoxs  that,  lor  xertieallx  polarized  radiation, 
the  scale  of  this  \  ariation  is  smalh'st  at  nc'ur-naclir 
angles  where  look  angle  is  almost  peqienclictila.'' 
to  the  surface  (Fig.  8)  (Kong  c't  al..  1979;  Fung  and 
Cihen.  1981;  Stogrxn,  1981;  Eppler  et  al..  1984). 

KRM.S  scans  (image  rows)  wc'cc  sumpic'cl  to 
c'Xtraet  the  middle  128  pixc-ls  (pixels  192-320)  so 
that  scene  nadir  falls  between  the  middle'  two 
columns  ofe'ach  te'sf  image  (columns  64  and  65). 
Look  angle'  at  the  extrc'ines  of  these  test  images 
is  approximately  13°  KRMS  data  analyzed  prc'- 
xiously  suggest  that  the  radiornc'tric  effc'ct  of 
change's  in  look  angle'  across  these  central  pixeds 
is  h'ss  than  2  K  for  first-xe-ar  sea  ice  (Fig.  Sb). 
Statistics  that  summarize  minimum,  maximum, 
and  mc'an  radiance  for  pixe'ls  across  the  tc'st  im 
agC'S  (Table  1)  corroborate  these  prexious  data 
and  sngge-sl  that  bias  across  the'  te'st  image's  is 
not  significant.  If  data  from  largc'r  atigh'S  xvere 
inclncle'd,  tlm  xariancc-  structure  displayed  at  the 
e'licls  of  eact>  roxx  could  be  significantly  ciifferc'nt 
from  that  in  the-  ce-ntral  peirtion  of  the  scan  line. 
This  xxould  xiolate-  .-Assumption  1  aboxe-,  xvhich 
recpiires  th  -t  the-  data  be  stationary  in  a  statistical 
sense. 

Scan-to-scan  change's  in  racliatice'  that  arise 
from  dilfe-rence-s  be-txxe-e'n  response-  characte-ristics 
of  the-  three  ante-nnas  xvere-  e-liminated  by  using 
data  de-rivecl  from  one-  of  th  three  KR.MS  antc'n- 
nas.  In  e-fle-c  t,  c  xe-ry  third  line-  of  the-  full  scene 
was  sarnple-d.  Tliis  re'sults  in  fore-shortc-ning  of 
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Figure  S  X'ertiealK  polari/.eel  hrightiu'ss  ter-mpt-ratun-  as  a 
liim-tioii  ol  inc-kle-iKe  aiiule-  .\)  as  pre-ciic  te'd  for  sea  wate-r. 
fresh  xvate-r  ice-,  aiie)  <lr\  snow  1)\  a  raeliome-tric;  inode! 
(from  Stours  n,  19S1)  and  B)  as  me-asiirt-d  from  KRMS  im- 
ai;es  of  first-sear  sea  iee-  (from  Kppler  et  al  .  19S4).  The 
sertieal  dotle-d  line-  in  .-X)  marks  the  imaging  limit  for  fnll- 
widtli  KRM.S  scenes.  The  image  limit  of  test  scenes  used 
he-re  is  approximafe-Ix  13°.  Tn  is  ambient  air  te-mperature 
and  7,.,  is  sk\-  temperature-  nse'd  in  the-  model  calculation. 


image  gc'ometry  by  a  factor  of  0.33,  and  reduction 
')f  overlap  between  adjacent  pixels  in  the  doxvn- 
column  direction  from  80%  to  40%.  Ovc'rlap. 
although  reducc'cl,  remains  sufficient  to  ensure 
statistical  dejjentlence  betwec'ti  adjacc'nt  pixels 
(Assumption  4  above).  Scene  geometrx'  and  overlap 
in  the  across-reivv  direction  rein,  ins  unaffected. 
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APPUCATION  OF  THE  WIENER  FILTER 

In  applying  the  Wiener  technique  to  the  three 
test  images  in  Figure  2,  we  evaluated  the  ett'ect 
of  altering  values  assigned  to  four  key  variables 
used  to  design  any  given  Wiener  filter:  1)  lag  (r 
in  Eqs.  (1)  and  (2)]  or  the  number  of  pixels  the 
original  scene  is  offset  to  define  the  desired  signal; 
2)  length  of  the  final  filter  \p  in  Eq.  (3)];  3)  number 
of  lines  (channels)  to  include  in  the  filter  (n, 
an  arbitrary  choice);  and  4)  degree  of  weighing 
applied  to  the  empirical  correlation  functions 
(both  autocorrelation  and  cross  correlation)  in  the 
estimate  of  the  true  correlation  (Robinson,  1983). 
In  the  present  situation,  another  option  is  to  de¬ 
fine  filters  that  are  applied  either  parallel  to  the 
scanning  direction  (across  the  image),  or  perpen¬ 
dicular  to  the  scan  direction  (down  the  image). 

The  evaluation  scheme  we  will  apply  involves 
mainly  visual  appraisal  of  the  processed  image. 
That  is,  does  one  combination  of  variables  "look” 
any  better  in  terms  of  presenting  desired  informa¬ 
tion  than  another  combination?  Statistics  such  as 
root  mean  square  are  calculated  that  are  useful 
for  relative  comparisons  (Robinson,  1983),  but  it 
is  ultimately  the  visual  clarification  of  the  textural 
components  in  each  scene  which  drives  this  analy¬ 
sis.  The  results  presented  below  were  selected 
from  many  combinations  of  lag,  length,  weight, 
and  number  of  included  lines  that  were  tried. 


Lag 

Figure  9  presents  the  original  three  scenes  and 
the  results  of  applying  a  filters  with  lags  of  1,  2, 
3,  and  5  pixels.  The  filter  length  was  held  at  a 
constant  length  of  25  pixels  for  all  scenes.  Addi¬ 
tionally,  three  channels  were  incorporated  in  the 
filter  design,  no  weighing  was  used,  and  the  im¬ 
ages  were  processed  parallel  to  the  scan  lines. 
These  images,  and  all  subsequent  images,  are 
presented  in  histogram  equalized  format  and  in¬ 
tensity  stretched  to  cover  the  pos.sibie  range  of 
intensities.  This  processing  enhances  subtle  fea¬ 
tures.  Figure  10  presents  the  original  images  mi¬ 
nus  the  processed  image,  also  equalized  and 
stretched  as  previously  described.  Such  presenta¬ 
tion  (Fig.  10)  is  useful  for  determining  what  effect, 
if  any,  a  particular  filter  had  on  the  original  image. 

Figures  9  and  10  show  that  increasing  the  lag 


between  the  original  and  desired  scene  removes 
apparent  noise  of  increasing  spatial  scale.  That  is, 
the  greater  the  lag,  the  greater  the  scale  of  noise 
that  is  minimized,  .additionally,  at  low  values  of 
lag  energy  associated  with  sharj)  intensify  changes 
is  smeared  over  a  larger  region  than  for  large  lags 
as  is  evident  in  the  first  column  of  images  that 
show  a  bright  feature  near  the  image  center. 
Figure  10  more  clearly  shows  the  enhancement 
of  progressively  longer  events  at  the  expense  of 
the  higher  frequency  information. 


Figure  1 1  presents  the  results  of  applying  filters 
with  lengths  of  3,  5,  7,  and  10  pixels.  The  lag  was 
held  to  1,  the  number  of  lines  (channels)  was  3, 
and  no  weights  were  applied.  Figure  12  displays 
the  original  minus  filtered  images  for  compari¬ 
son.  Visually,  the  advantage  of  using  longer  than 
3-pi.xel-length  filters  is  difficult  to  ascertain  in  this 
case  as  it  appears  that  the  majority  of  the  noise 
reduction  w^as  accomplished  with  shorter  filters. 
It  is  interesting  to  note  that  filters  in  the  range  of 
3-5  pi.xels  correspond  roughly  to  the  amount  of 
average  signal  overlap  contained  in  each  scan, 
which  is  in  alignment  with  the  assumptions  associ¬ 
ated  with  the  application  of  Wiener  theory.  Very 
little  improvement  is  realized  with  increased  filter 
length,  suggesting  that,  at  least  on  the  scale  of 
variability  shown  in  these  images,  significant  en¬ 
hancement  can  be  obtained  using  relatively  short 
filters.  However,  examination  of  the  normalized 
rov't  mean  square  error  associated  with  a  given 
filter  length  (Table  3)  suggests  that,  at  least  rela¬ 
tive  to  the  case  described  above,  the  majority  of 
the  decrease  in  this  error  is  accomplished  with  a 
filter  length  of  12  pixels  and  under.  This  is  more 
clearly  demonstrated  in  Table  4,  which  depicts 
the  amount  of  change  in  the  normalized  root- 
mean-square  error  as  filter  length  is  increased. 
Further  examination  of  Tables  3  and  4  indicate 
that  the  same  is  generally  tnie  for  lags  of  2  and 
3  pixels  whereas  filters  defined  for  longer  lags 
would  benefit  by  an  increased  filter  length. 


Changing  the  number  of  lines  (channels)  used  to 
define  the  Wiener  filter  produced  results  dis¬ 
played  in  Figures  13  and  14.  Th('  number  of  lines 


Length 


Lines  (Channels) 


14  Full  and  Eppler 


ABC 


Original 


Lag  1 


Lag  2 


Lag  3 


Lag  5 


Figure  9.  Effects  of  using  lags  of  different  magnitude  in  Wiener-filtered  images.  Uiifilferec) 
images  are  shown  in  the  first  row  (original}.  Images  filtered  with  shifts  of  1,  2.  3.  and  .5 
pixels  in  the  horizontal  direction  {lag=  1, 2,3.5)  are  shown  in  successive  rows,  filter  length 
was  held  constant  at  25  pixels,  three  adjacent  lines  were  used  to  define  a  filter  operator, 
no  weighing  applied,  and  processing  was  parallel  to  the  scan  direction  (horizontal  in  these 
images).  Note  that  as  lag  increases,  longer  wavelength  feattires  are  enhanced  while  shorter 
wavelength  features  are  minimized. 
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Figure  10.  Difference  images  derived  by  subtracting  the  filtered  images  in  Figure  9  from 
the  unfiltered  images  (first  row).  Features  that  appear  in  these  different  images  represent 
information  removed  from  the  unfiltered  scenes. 


used  in  each  row  of  images  varied  from  1  to  10. 
The  lag  was  1  for  all  images,  the  length  of  the 
filter  was  25,  and  no  weights  were  applied.  Since 
Wiener  theory  relies  on  the  autocorrelation  and 
cross-correlation  functions,  highly  correlatable 
events  will  tend  to  dominate  design  of  the  opti¬ 


mum  filter.  Such  ca.ses  arise  in  our  data  along  the 
cage  of  clearly  defined  features,  in  the  vicinity  of 
a  textural  change,  or  at  abrupt  changes  in  inten¬ 
sity.  The  net  result  is  that  the  greater  the  number 
of  lines  used  in  the  analysis,  the  blockier  the 
filtered  image  appears  to  be  in  Figure  14.  Also 
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Length  =  5 


Length  =  7 


Length  =  10 


Figure  11.  Effect  of  using  Wiener  filters  of  different  length  (3,  5,  7,  and  10  pixels).  Lag 
was  one  pixel,  three  lines  were  used,  and  no  weighing  function  was  incorporated.  Note 
the  lack  of  visually  significant  changes  as  increasingly  longer  filter  operators  are  applied. 


noted  in  these  images  is  the  tendency  of  filters 
composed  of  larger  numbers  of  lines  to  enhance 
the  more  dominant,  correlatable  events  in  a  scene, 
which  are  defined  by  a  larger  number  of  lines. 
This  results  in  less  enhancement  of  smaller-scale 
features  which  are  defined  by  relatively  few  lines. 
Using  two  or  three  channels  in  the  filter  design 
would  roughly  correspond  to  the  average  amount 


of  overlap  among  adjacent  scans,  which,  as  noted 
before,  corresponds  to  the  assumptions  which 
underlay  this  analysis. 

Weight 

Changing  the  weighing  function  (Figs.  15  and  16) 
produced  no  visually  significant  differences.  The 
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Figure  12.  Diflference  imagec  derived  from  subtracting  the  filtered  images  in  Figure  1 1 
from  the  unfiltered  images  (first  row).  Only  subtle  changes  are  observed  as  filter  length  in- 
creases. 


lag  was  1,  the  length  was  25,  and  three  lines 
were  included  in  the  filter.  The  weighing  function 
values  displayed  in  these  images  include  0  (no 
weights),  10  (maximum  weighing),  and  30  and  50 
(intermediate  weights)  [weighing  scale  defined  in 
Robinson  (1983)].  The  conclusion  we  draw  from 
these  images  is  that  weighing  for  the  correlation 
function  is  not  a  visually  significant  variable  to 


consider  for  this  type  of  image  data.  Further  analy¬ 
sis  may  be  warranted  when  the  image  displays 
different  variance-covariance  structures. 

Scan  Direction 

Finally,  scenes  were  processed  fjoth  parallel  to 
the  scan  direction  and  perpendicular  to  the  scan 


18  Full  and  Eppler 


Table  3.  Filter  Length  and  Norinidized 
Root-Mean-Square  Error  Values  for  Syste  ms  of  1,  2, 

3,  4,  and  5  Lag  Pixels  Presented  for  Filter  Lengths  up 
to  25  Pixels” 


Normalized  Root  Mean  Square 


Ijengfh 

Lag=l 

II 

Lag  =  3 

Lag  =  4 

Ug  =  5 

1 

0.339 

0.372 

0.425 

0.492 

0,579 

2 

0.304 

0.343 

0.408 

0.476 

0.568 

3 

0,293 

0.337 

0.393 

0.470 

0,559 

4 

0.288 

0.319 

0.386 

0.455 

0..539 

■5 

0.269 

0,311 

0.37.' 

0.436 

0,,504 

6 

0.266 

0.305 

0.362 

0.419 

0474 

t 

0.260 

0.295 

0.345 

0.388 

0,458 

8 

0.233 

0.283 

0.321 

0.374 

0.443 

9 

0.241 

0.261 

0.303 

0,357 

0.422 

10 

0.230 

0.249 

0.298 

0.349 

0.404 

11 

0.218 

0,244 

0.290 

0.336 

0.394 

12 

0.210 

0.240 

0.280 

0.332 

0..389 

13 

0.208 

0,233 

0,273 

0.325 

0.383 

14 

0.202 

0.227 

0,268 

0.320 

0.375 

15 

0.197 

0.224 

0.262 

0.313 

0.363 

16 

0.193 

0.219 

0.256 

0..300 

0.350 

17 

0.188 

0.215 

0.250 

0.294 

0.342 

18 

0.184 

0.210 

0.244 

0.282 

0.,330 

19 

0.178 

0.205 

0.237 

0.273 

0.316 

20 

0.176 

0.199 

0.227 

0.262 

0.297 

21 

0,171 

0.192 

0.218 

0.244 

0,291 

22 

0.167 

0.185 

0.207 

0.234 

0.270 

23 

0.161 

0.178 

0.198 

0.226 

0,260 

24 

0,153 

0,171 

0.191 

0.214 

0.238 

2.5 

0.148 

0,167 

0.186 

0.201 

0.226 

"  These  values  are  recorded  for  the  first  set  of  scan  lines  in  the  first 
image  and  are  typical  for  the  majority  of  the  scans  in  each  image. 
There  were  three  scan  lines  used  to  define  each  filter. 


Table  4.  Change  in  \orrnali/,e<i  Root-.Mean-Squure 
Error  Values  for  Filters  of  Increasing  l-ength  in  .Systems  of 
I.  2.  3.  4,  and  5  Lag  Pixels” 


Step 

Cdiange  in  Nonnalized  Root  Mean  Sipiare 

II 

lMg^2 

li 

••sa 

Lag  =  4 

iMg^S 

1-2 

0.035 

0.029 

0.017 

0,016 

0  01 1 

2-3 

0.011 

01)06 

0  015 

0  (  M  )6 

0009 

3-4 

0.(K)5 

0  018 

0.007 

0.015 

0  020 

4-5 

0.019 

0,008 

0.013 

0.019 

0.035 

5-f> 

0,003 

0006 

0  011 

OOIT 

0  030 

6-7 

0006 

0.010 

0.017 

0.031 

0.0  Ih 

7-8 

0.(M)7 

0.012 

0,024 

0  014 

0.015 

8-9 

0,012 

0.022 

0,018 

0017 

0.021 

9-10 

0-01 1 

0.012 

0.005 

0.008 

0  018 

10-11 

0  012 

0.005 

0  006 

0,013 

OOlO 

11-12 

0.008 

0,(M)4 

0,010 

0.(«)4 

0  005 

12-13 

0002 

0007 

0  007 

0  007 

0006 

13-14 

0.006 

0.006 

0.005 

0(8)5 

0,(K)8 

14-15 

0.005 

0.003 

0008 

0,007 

0,012 

1.5-16 

0.(M)4 

0005 

O.OW) 

0  013 

0013 

16-17 

0.005 

0.004 

0.006 

0(8)6 

0(K)8 

17-18 

0.004 

0.00,5 

0006 

0.012 

0012 

18-19 

0,006 

0.005 

0.007 

0.009 

0.014 

19-20 

0,002 

0.006 

0  010 

0  01 1 

0.019 

20-21 

0,005 

0.007 

0009 

0,018 

0,006 

21-22 

0,004 

0.007 

0,011 

0.010 

0,021 

22-23 

0,006 

0.007 

0,009 

{),(M)8 

0010 

23-24 

0.008 

0.007 

0.007 

0,012 

0.022 

24-25 

0.005 

0.004 

0.005 

0.013 

0  012 

“  These  values  were  ealeiilated  from  Table  : 

3.  The  step 

values  refer 

to  the  fact  that  one  ro«  of  Table  .3  repr<*senting  a  given  filter  length 
is  subtracted  from  the  previous  row,  whii  h  reflects  a  filter  length  I 
pixel  shorter.  Larger  values  in  the  change  columns  iudicafe  relatisely 
larger  changes  in  root-inean-srpiare  salues  with  increasing  filter 
length. 


direction.  No  visually  significant  differences  be¬ 
tween  images  were  observed  except  that  the 
scenes  processed  parallel  to  the  scan  lines  seemed 
a  little  sharper  than  those  processed  perpendicu¬ 
lar  to  the  scan  lines.  This  may  be  due  to  the  fact 
that  there  is  a  larger  correlation  between  adjacent 
pixels  along  a  scan  than  there  is  between  adjacent 
pixels  in  two  scans  simply  based  on  amount  of 
oversampling  by  the  instrument. 

DISCUSSION 

The  application  of  multichannel  Wiener  theory 
to  passive  microwave  images  of  first-year  sea  ice 
has  produced  visually  sharper  images.  That  is,  the 
filtered  images  appear  to  be  more  focused  than 
the  original  images.  The  advantages  of  the  Wie¬ 
ner-processed  images  are  that  textural  features 
are  enhanced,  and  high  frequency  noise  is  re¬ 
duced.  The  disadvantages  are  that  this  filter  in¬ 


duces  a  slight  blocky  grain  to  the  image  and  can 
produce  a  type  of  ringing  along  scan  lines  which 
traverse  sharp  intensity  contrasts  (discontinuous 
variability  in  the  scans).  The  usual  practice  of 
applying  tapered  filter  weights  (Rush,  1990)  to 
remove  this  ringing  is  not  normally  performed 
with  the  Wiener  filter.  Rather,  the  existence  of 
such  ringing  is  probably  indicative  of  a  local  pixel 
region  whose  statistical  structure  differs  markedly 
from  the  regional  structure.  Defining  regions  with 
similar  statistical  structure  before  the  application 
of  Wiener  theory  might  allow  the  definition  of  a 
single  filter  that  can  be  used  for  this  region.  This 
filter  might  be  defined  from  previous  experience 
or  simply  be  the  product  of  the  first  set  of  scan 
lines  encountered  in  the  defined  region.  In  our 
case,  the  advantages  outweigh  the  disadvantages 
because  ice  textural  features  that  occur  on  inter¬ 
mediate  scales  are  enhanced  which  is  the  primary 
objective  of  this  study.  A  reasonable  filter  to  apply 
in  our  case  would  consist  of  two  or  three  channels, 
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Figure  13.  Effect  of  using  Wiener  filters  that  are  composed  of  different  numbers  of  scan 
lines  (2,  3,  4,  and  10).  In  each  case,  lag  was  one,  filter  length  was  25,  and  no  weighing 
function  was  used.  Image  texture  becomes  increasingly  blocky  as  the  numljer  of  lines  in- 
creases. 


have  a  lag  of  1  pixel,  a  length  of  approximately 
12-14  pixels,  and  no  correlation  weighing.  This 
image  would  be  filtered  parallel  to  the  scan  direc¬ 
tion.  This  combination  is  based  on  the  minimum 
filter  needed  to  show  significant  visual  improve¬ 
ment.  More  complex  filters  have  not  been  shown 
to  increase  visual  enhancement  significantly. 


In  terms  of  practical  application,  however, 
Wiener  theory  does  have  some  drawbacks.  For 
example,  the  Wiener  algoritlun,  as  implemented 
by  us  in  this  study,  defines  a  new  set  of  filters  for 
each  collection  of  scan  lines.  In  a  128  x  128  scene, 
a  three-channel  filter  would  involve  129  separate 
applications  of  the  algorithm.  This  approach,  while 
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Figure  14.  Difference  images  derived  by  subtracting  the  filtered  images  in  Figure  13  from 
the  unfiltered  images  (first  row). 


•V 
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Figure  15.  Effect  of  weighing  on  Wiener-filtered  images:  no  weighing  (0),  maximum 
weighing  (10),  and  intermediate  weights  (30  and  50).  Lag  was  one  pixel,  three  lines  used 
for  filter  definition,  and  filter  length  was  25  pixels  in  all  cases.  Vistiallj’  significant  changes 
are  not  observed  as  the  weighing  function  varies. 
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Figure  16.  Difference  images  derived  by  subtracting  the  filtered  images  in  Figure  15  from 
the  unfiltered  images  (first  row).  These  difference  images  generally  support  the  conclusion 
that  no  visually  significant  differences  appear  when  various  weights  are  applied  in  the  fil¬ 
ter  design. 


computationally  cumbersome,  was  chosen  in  or¬ 
der  to  allow  the  filtering  scheme  to  be  more  sen¬ 
sitive  to  changes  in  altitude,  plane  position,  and 
changes  in  ice  type.  In  such  an  approach  a  single 
filter  appropriate  for  a  specific  ice-type  or  ice 
texture  would  be  defined.  This  would  result  in 
significant  computational  advantages.  Such  filters 
could  not,  however,  be  applied  across  abrupt 


changes  in  intensity  such  as  those  associated  with 
different  ice  types.  The  result  of  such  an  applica¬ 
tion  would  be  the  introduction  of  a  ringing  in  the 
region  of  the  abrupt  change.  If  such  areas  were 
included  in  images  used  to  design  the  filter,  the 
resultant  filter  would  enhance  the  largest  and 
most  correlatable  events  at  the  expense  of  smaller 
features,  such  as  the  more  subtle  texture  observed 
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in  first-year  ice.  This  is  a  CGinn.on  problem  associ¬ 
ated  with  most  variance  algorithms  and  has  re¬ 
ceived  extensive  discussion  in  other  fields  such  as 
seismic  analysis  (Wadsworth  et  ah,  1953). 

Passive  microtcave  imagery  used  in  this  study  was  acquired 
in  conjunction  with  the  Joint  Navy-S'ASA  SSM/I  Sea  Ice 
Validation  Experiment  usmg  funds  provided  both  by  the 
Oceanographer  of  the  Navy  (OP-096,  project  element  63704N) 
through  the  Satellite  Applications  Technology  (SAT)  Program, 
A.E.  Pressman  program  manager,  and  by  .V.4SA  through  the 
Snow  and  Ice  Processes  Program,  R.  Thomas  program  man¬ 
ager.  Data  analysis  teas  funded  under  Department  of  Defense 
program  element  0601153N  thnmgh  the  Sea  Ice  Classification 
Project  and  through  a  grant  from  these  funds  to  W.  E.  E. 
M'  E.  F.  also  received  support  under  a  Summer  Faculty  Re¬ 
search  Fellowship  awarded  by  the  U.S.  Navtj-American  Soci¬ 
ety  for  Engineering  Education  Fellowship  Program  through 
NOARL's  Basic  Research  Management  Office  (Code  IN),  U.  E. 
Morris  head.  This  work  represents  NOARL  Contribution  No. 
JA:332:()19:92. 
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